Engineering light localization in a fractal waveguide network 
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We present an exact analytical method of engineering the localization of electromagnetic waves in 
a fractal waveguide network. It is shown that, a countable infinity of localized electromagnetic modes 
with a multitude of localization lengths can exist in a Vicsek fractal geometry built with diamond 
shaped monomode waveguides as the 'unit cells'. The family of localized modes form clusters of 
increasing size. The length scale at which the onset of localization for each mode takes place can 
be engineered at will, following a well defined prescription developed within the framework of a real 
space renormalization group. The scheme leads to an exact evaluation of the wave vector for every 
such localized state, a task that is non-trivial, if not impossible for any random or deterministically 
disordered waveguide network. 

PACS numbers: 42.25.Dd, 42.82.Et, 72.15.Rn 



I. INTRODUCTION 

Since its inception, the phenomenon of Anderson lo- 
calization in a disordered lattice [l| has been found to be 
ubiquitous in diverse fields of condensed matter physics 
and materials science. While the popular domain of inter- 
est is related to electronic systems, where the quantum 
interference plays a pivotal role in localizing electronic 
eigenstates in presence of disorder , the effect is by 
no means, confined to electrons, and extend over a vari- 
ety of phenomena ranging from spin freezing in one di- 
mensional semiconductors @, localization in optical lat- 
tices [Tol - IT3 ] . or the localization of matter waves (cold 
atoms forming Bose-Einstein condensates) [I]! [l4j, to 
name a few. Incidentally, the latter has recently been ob- 
served experimentally in one dimensional matter waveg- 
uides, where the random potential is generated by laser 
speckles (l5j . 

In 1984 Sajeev John pointed out that the idea of lo- 
calization goes far beyond the electronic systems, and 
is actually a general phenomenon common to any wave 
propagation in systems with disorder jTfj| . Anderson fol- 
lowed with a seminal paper considering the idea of lo- 
calization of classical waves, in an attempt to work out 
the theory of white paints [l7|- The field gathered mo- 
mentum in the last couple of decades and a considerable 
volume of literature related to the localization of classi- 
cal waves, particularly light, ultrasound and microwave 
is now in existence [18l426l |. 

The study of localization of light in disordered me- 
dia has been patronized by the discovery of the photonic 
band gap (PBG) materials [27i - [29| . These systems ex- 
hibit gaps in the frequency spectrum in which propa- 
gation of electromagnetic waves is forbidden. This has 
important implications in both fundamental science and 
technological applications. 

Photonic gaps, apart from materials with a large di- 
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electric constant, can also be observed in waveguide 
networks, as proposed by several groups over the past 
years |3fj| - l38l |. Anderson localized eigenmodes are ob- 
served inside the photonic gaps and excellent agreement 
between theory and experiments has been obtained [3(| ■ 
The network models are able to localize a propagating 
wave by virtue of the geometrical arrangement of the 
waveguide segments. Of particular interest are a wide 
variety of models based on waveguide networks designed 
following a deterministic fractal geometry [34U381 ] , where 
the gaps result from the typical topology exhibited the hi- 
erarchical arrangement of the waveguide segments. The 
present communication also deals with a hierarchically 
designed (fractal geometry) waveguide network, but ad- 
dresses a deeper fundamental question regarding wave 
localization in such systems, as explained below. 

Fractals or hierarchical geometries in general cause 
an excitation to localize [39j-|42j]. The energy spectrum 
turns out to be singular continuous (_{__ . with a gap in 
the vicinity of every energy. Nevertheless, there can 
be a countable infinity of extended eigenfunctions with 
high (or even, perfect) transmittivity, even though there 
is no translational order in such systems. Once again, 
this is true for electrons (43l - |49l |. and classical waves as 
well [35|,[36|. A curious point, apparently gone unnoticed 
or un-appreciated so far is that, while a precise determi- 
nation of the eigenvalues corresponding to the extended 
eigenmodes is possible in the above cases of hierarchi- 
cally grown fractal networks, the task seems to be prac- 
tically impossible when it comes to an exact evaluation 
of eigenvalues of the localized modes in such hierarchical 
systems in their thermodynamic limits. Direct diagonal- 
ization of the Hamiltonian (in the electronic case) , or an 
exact numerical solution of the wave equation doesn't 
help, as the overall character of the spectrum in all the 
cases is highly fragmented, and the eigenvalues obtained 
from a finite sized network are likely to slip away from 
the spectrum once we go over to a higher generation. 
This problem has recently been addressed in the context 
of electron localization in fractal space I50|. and we carry 
forward the central idea floated in Ref. |5C| to evaluate 
the exact wavelengths (wave vectors) of electromagnetic 
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waves that can be localized at will in a properly designed 
hierarchical single channel network. To the best of our 
knowledge, this issue remains unaddrcssed so far in the 
literature in the context of localization of wave, light in 
particular. 

We design a Vicsek fractal network [5l| consisting of 
diamond shaped loops each arm of which mimics a single 
mode linear waveguide [see Fig.[T]. While examination of 
the localized mode eigenvalues and the nature of localiza- 
tion are indeed the major factors motivating this work, 
other interests in such a study are related to the general 
spectral character and classical wave transport in these 
systems. As one can easily appreciate, the geometry of 
a diamond- Vicsek network provides an interesting culmi- 
nation of the 'open' character of a typical Vicsek pattern 
and closed loops at shorter scales of length. This is in 
marked contrast to the much studied Sierpinski gasket 
waveguide network [34|, which is a closed structure, or 
to the other deterministic waveguide networks [36|, [3?| ■ 
The presence of the loops generates a possibility of an 
effectively long ranged propagation of waves between the 
various vertices, and its effect on the localization or de- 
localization of waves is worth studying. 

We find interesting results. For an infinite hierarchical 
geometry, such as presented above, a countable infinity of 
cigcnmodcs with a multitude of localization lengths can 
be precisely detected. One can work out an exact mathe- 
matical prescription to specify the length scale at which 
the onset of localization begins. The localization can in 
principle, be delayed (staggered) in position space and the 
corresponding wave vectors (or, wavelengths) can be ex- 
actly evaluated following the same prescription based on 
a real space renormalization group (RSRG) method (52j. 
In addition, it is shown that for a given set of parame- 
ters, the center of the spectrum corresponds to a perfectly 
extended eigenmode, with the parameters describing the 
system exhibiting a fixed point behavior. This central ex- 
tended mode is flanked on either side by localized wave 
functions with a hierarchy of localization lengths. 

In what follows we describe the results. In Sec. [Til 
the model and the mathematical method of handling the 
problem are presented. Section [TTT] discusses the results 
and their analyses, and in Sec. IIVI we draw our conclu- 
sions . 



II. THE MODEL AND THE METHOD 

A. The wave equation and its discretization 

We have considered a waveguide network formed by 
waveguide segments having the same lengths arranged 
in a Vicsek fractal geometry Each segment has a 

single channel for wave propagation. The wave function 
tpij between any two nodal points i and j satisfies the 
wave equation: 



(1) 

where uj is the frequency of the wave, c is the speed of 
wave propagation and x is the distance measured from 



the ith node. The above equation has the solution of the 
form [H, 



4>tj{x) = fa 



sm[k(£ij — x)] 
sm(k£ij) 



sin(for) 
sin(k£ij) 



(2) 



where k — u)/c, £ij is the length of the segment between 
the node i and j, and ipi and tpj are the values of the 
wave function at the ith and jth nodes respectively. The 
flux conservation condition: 



E 







(3) 



x=0 



where the summation j is over all the nodes linked di- 
rectly to i, leads to a discretized version of Eq. ([!]) f54j . 
viz., 



ipi ^2 cot % + E sin 







(4) 



wher 



k£t 



ka, a being the constant length of a 



waveguide, as considered in all the calculations which fol- 
low. Eq. Q resembles a tight binding difference equation 
depicting the propagation of non-interacting electrons in 
a lattice, viz., 



(e - ei)ipi = y] Ujipj 



(5) 



with the 'electron energy' E can be put equal to 2 cos ka, 
the 'on-site potential' q = 2 cos ka + c °t % an d 'hop- 
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ping matrix element' Uj = 1/ shifty. It should be ap- 
preciated that the choice of E = 2 cos ka is completely 
arbitrary, and this does not affect the final results in any 
way. As seen in Fig. [TJ the Vicsek waveguide network has 
two types of nodal points, viz., type A (having two neigh- 
boring nodal points) and type B (having four neighboring 
nodal points). Accordingly, the equivalent 'on-site poten- 
tials' are assigned two values, viz., = 2 cos ka+2 cot ka 
and cb = 2 cos ka + 4 cot ka respectively. The 'overlap 
integral' along an arm of the waveguide is t = 1/sinfca. 
There is no second neighbor tunneling of the wave to 
begin with. But it will grow with RSRG. 

We now proceed to describe the physics of wave propa- 
gation in such a fractal geometry by exploiting this exact 
analogy with the corresponding electronic problem. 



B. The RSRG scheme 

A renormalized version of the fractal network is easily 
obtained by decimating a subset of vertices from the orig- 
inal geometry. This implies one has to eliminate a subset 
of the wave amplitudes from the difference equations ([5]) 
in terms of the surviving vertices [see Fig. G^b)]. This 
results in the following set of recursion relations for the 
system parameters. 



FIG. 1: (Color online) (a) Schematic view of the second generation of an infinite diamond- Vicsek waveguide network. The 
vertices at the junction of two waveguides are named as l A' in the text, while those in the bulk (at the crossing of four 
waveguides) are termed 'B\ (b) Renormalized version of (a) with the dotted lines indicating the diagonal 'hopping' which is 
generated due to renormalization. 
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where ^ = 2t 2 (E-E A -X)/A with A = (E — e A — X)(E — 
e b — X) — 2t 2 , and A = at the beginning. 

The scaling generates an effective second neighbor 
'hopping' (overlap) A', as is obvious from the above set of 
recursion relations and the diagramQJb). It is important 
to appreciate that, the growth of a A', unlike the elec- 
tronic case does not have any physical significance 
here, as there is no real 'overlap' of the wave amplitudes. 
The light is essentially confined within the waveguide. 

These recursion relations will now be used to obtain in- 
formation about the local density of eigenmodes at spe- 



cific sites of the system, and the localized or extended 
character of the modes, as discussed below. 



III. RESULTS AND DISCUSSION 

A. Local density of eigenmodes 

As already stated, the exact mapping of the wave equa- 
tion on to a discrete Schrodinger type equation allows us 
to extract information about the density of electromag- 
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netic modes through a Green's function analysis [52| ■ We 
present in Fig. [2] the density of modes at a B- vertex, 
which is given by, 



(ka) 



lim 



lm{G {B) (ka + ii])} 



(7) 



The distribution of eigenmodes, plotted within ka < < 
27r, shows clusters of non-zero values over a finite range of 
the wave vector fc, and is found to be symmetric around 
ka = tt/2. The fragmented Cantor- like character, typi- 
cal signature of a fractal spectrum, is apparent. An idea 
about the character of the eigenmodes can be obtained 
by observing the flow of t under successive RSRG itera- 
tions. In general, for an arbitrary value of ka, for which 
the density of modes is non-zero, — > as the number 
of iterations n increases, implyin g a localized character 
of the corresponding eigenmode |43t l52j . However, for 
ka = (2m + l)%/2, both the first and the second neigh- 
bor 'hopping integrals' t and A remain non-zero for an 
indefinite number of iterations. In fact, we observe a 
one cycle fixed point of the entire parameter space, viz., 
{e' A , e' B , t', A'} = {eA,£B,t, A}. The fact that t and A re- 
main finite at all stages of RSRG implies that there is 
a non-zero overlap between the wave amplitudes at all 
scales of length, and the corresponding mode is an ex- 
tended one. 

The neighborhood of ka = (2m + l)7r/2 has been 
scanned minutely. The self similarity of the spectrum is 
always seen with dense patches of eigenvalues clustered 
throughout the intervals. For many of these eigenvalues 
t and A remain finite under successive decimation for a 
large number of steps. This indicates that every local 
band center at ka = (2m + l)7r/2 is flanked either by 
extended modes or at least, eigenmodes with very large 
localization lengths. 
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FIG. 2: (Color online) Local density of eigenmodes at a 
bulk B-type vertex of an infinite diamond- Vicsek waveguide 
network. 



B. Explicit construction of localized modes 

The deterministic Vicsek fractal waveguide is self- 
similar at all scales of length. This feature allows us to ex- 
plicitly construct a special distribution of wave functions 
by suitably exploiting the difference equations Eq. ([5]). 
These 'special' eigenmodes are localized and extend over 
clusters of single channel waveguide segments of various 
sizes. The planar extent of such clusters depends on the 
eigenvalue corresponding to the localized mode, and can 
be small or enormous, depending on the wavelength (or, 
wave vector). The construction is similar to our recent 
work [50| on electronic states. 

To elaborate, let us set 



E = e B (n) - 2X(n) 



(8) 



where, n refers to the stage of renormalization. This is 
in general, a polynomial equation in E (and hence, in k). 
The zeros of this equation will be the allowed wave vec- 
tors for the infinite system if, and only if, with them, one 
can satisfy Eq. locally at every vertex of the network. 
This task can be accomplished by trying to draw a non- 
trivial distribution of amplitudes for a value of E (ka) 
obtained from Eq. (|5J| on the undecimated vertices of an 
n-step renormalized network, and then trying to figure 
out the amplitude distribution on the original waveguide 
structure at the bare length scale. This can indeed be 
done, as we demonstrate in Fig. G^a) for n = 1. For 
n = 1 Eq. © reduces to 



cos ka(2 cos 2ka — I) = 



(9) 



Roots of the above equation are ka = ±7r/2 and ka = 
±7r/6. ka = n/2 (or, equivalently, —n/2) is of course, 
the extended mode. The root ka = ir/6 leads to the con- 
struction of wave amplitudes as shown in Fig. [3] It is not 
difficult to extend the construction depicted in Fig. G3a) 
even to a network of an arbitrarily large size, where the 
end vertices are not actually visible. We are still able 
to satisfy Eq. §5§ locally at every vertex while drawing 
this distribution and thus, k = n/6a definitely belongs to 
the spectrum of the infinite system, a fact that has been 
cross-checked by evaluation the LDOS at the A- and the 
B-sites at this special value of k. We get a stable, finite 
value of the LDOS which supports our argument above. 

In Fig. [SJa) we show the distribution of amplitudes on 
the central cluster of an infinite diamond- Vicsek hierar- 
chical network for ka = it/ 6. The deep red arms connect 
network vertices where the wave amplitudes are non-zero, 
and thus these arms are the brightest looking ones as far 
as the distribution of light intensity is concerned. The 
black lines represent waveguides which will appear com- 
pletely dark, as the wave amplitudes at their vertices will 
have to be zero in order to satisfy Eq. ([5]) . There will be 
arms connecting one vertex with zero amplitude and an- 
other with a non-zero one. These are depicted by a lighter 
shed of red, and will 'glow' with less intensity compared 
to the deep red ones. The distribution of intensity in any 
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FIG. 3: (Color online) (a) Distribution of amplitudes of the wave function for ka — n/6 (obtained by solving Eq. © for 
n — 1) on a second generation network. The dark shaded plaquettes with black lines at the boundary embrace network vertices 
with zero amplitude. The deep red waveguide segments (covering a dark red shaded region) connect vertices with non-zero 
amplitude, and 'glow' with maximum intensity. The lighter red lines represent waveguide segments which have an intermediate 
intensity profile. The amplitudes of the wave are marked by the numbers ±1, 0, ±a and ±/3 respectively, where a = — V3/2 
and P — 1/2. (b) The distribution of wave amplitudes for ka — tt/6 on a third generation network. The thick blue highlighted 
lines represent one- and two-step renormalized lattices in (a) and (b) respectively. 



arm (apart from the black ones) is by no means uniform. 
The colors just represent the fact that the intensity is 
non- zero. The significant observation is that, clusters of 
non-zero amplitude span over a finite distance, but ulti- 
mately get 'decoupled' from each other on a larger scale 
of length. This can be appreciated if we look at Fig.[3Jb) 
which is a larger version of the previous figure. The red 
shaded clusters are distributed along the principal X- 
and y-axes, but are separated from each other beyond 
a certain extent by light red boxes. The black clusters 
representing amplitude- voids are now seen to span larger 
spatial distances. A similar construction is possible for 
ka = — 7r/6 which is another solution of Eq. © for n = 1, 
but this does not have any additional significance. In 
terms of light, the entire hierarchical geometry will have 
an appearance where light will be localized with higher in- 
tensity at certain clusters of waveguides, decoupled from 
each other by completely dark patches. 

It is apparent from the above discussion that the eigen- 
function corresponding to ka = ±7r/6 will be localized in 
the fractal space. This is easily re-confirmcd by studying 
the evolution of the hopping integrals under successive 
RSRG steps. The hopping integrals t and A (zero ini- 
tially, but grows later) remain non-zero at the first stage 
of RSRG (that is, n = 1), indicating that the nearest 
neighboring sites on a one step renormalized lattice will 
have a non-zero overlap of the wave functions. They start 
decaying for n > 2 with the decay in A(n) taking place 



at a much slower rate compared to t(n). This indicates 
that over larger scale of length the corresponding states 
are localized, but the effect is a weak one. 

C. The staggering effect 

The previous observation immediately leads to an in- 
novative way of exactly determining the wave vectors 
(wavelengths) corresponding to localized wave functions 
on such a deterministic geometry at an arbitrary scale of 
length. We do it using the following method. 

We can solve Eq. © in principle, to get the desired de- 
values for any n. For example, we have done it explicitly 
for n — 2. The roots are obtained from the equation, 

cos ka(2 cos 2ka — l)(cos2fca — 2cos8fca) = (10) 

and are given in Table [I] In every case, on beginning 
the RSRG iteration with the wave vector k chosen arbi- 
trarily from the above set, the nearest neighbor overlap 
integral t, and the diagonal one, viz., A remain non zero 
at least up to that specific n-th stage of renormalization. 
After that, as the RSRG progress, the hoppings flow to 
zero with A dominating over t at every step of renormal- 
ization. This implies that for any such fc-value one can 
draw a non-trivial distribution of wave- amplitudes on the 
renormalized fractal network. When mapped back on to 
the original lattice the amplitudes will be found to span 
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n 


Values of ka 


1 


±1.570796, ±0.523599 


2 


±1.570796, ±0.523599, 
±2.617994, ±1.823731, 
±1.317862, ±2.132603, 
±1.008990, ±2.531027, 
±0.610566, ±3.008153, 
±0.133440 



by the well known formula [5f 



TABLE I: Values of ka obtained from the Eq. 
and n — 2 respectively. 



for n 



clusters of increasing size. The exact size of the spanning 
clusters will be determined by the value of n. The size, 
for example, with n = 2 exceeds that for n = 1. 

The spanning clusters finally get decoupled from sim- 
ilar clusters when one looks at the distribution over a 
large enough network. Speaking in terms of the red and 
grey shaded zones, it should be appreciated that the size 
of the red zone is much bigger for n = 2 compared to the 
n = 1 case. We refer the reader to Ref. [5jJ for under- 
standing the result. 

It is now obvious that higher the value of n, more 
will be the number of roots of the polynomial equation 
Eq. ([5]). The roots have a nice nesting property. The 
roots obtained from any (n — l)-th stage are found in- 
cluded in the solution for the n-th stage (see Table HJ). 
The additional roots obtained at the n-th level over the 
existing roots from the (n — 1)- th level keep the overlap 
integrals non-zero up to that specific n-th RSRG step. 
Beyond this step, the overlap finally starts to weaken 
in magnitude. This immediately implies that for larger 
values of n, the clusters of non- zero amplitudes span 
wider fractal space, and the localization effect begins 
much later. That is, on-set of localization can be delayed 
in space, at one's will, by choosing to solve Eq. (J5J for 
larger value of n. We can thus bring a staggering effect 
on the localization of light in such a fractal waveguide 
network. 



D. Transmission of electromagnetic waves 



T 



4 sin kc 



(11) 



where V x = [(M u - M 21 ) 
V 2 = [(Mix ±M 22 )sinfca] 



Mjj are given by, Mxx 



(E-U) 



Mai 



-M12, M 22 



(Mn — M 22 ) cosfca] and 
The matrix elements 

(E-u) 2 _ r 

TU i, ' 
ti 



M 12 
and cos ka 



(E — ei)/2ti, 'a' being the lattice constant and is taken 
to be equal to unity throughout the calculation. 
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FIG. 4: (Color online) Transmission characteristics for a 
3rd generation Vicsek waveguide network. The length of each 
monomode segment a is chosen as unity. 

In Fig. [51 we have shown the two-terminal transmis- 
sion characteristics for a 3rd generation waveguide net- 
work system. The transmission spectrum is, as expected, 
full of gaps, with resonant transmission exhibited at cer- 
tain fc-values. With increasing generation, the resonances 
become rarer and the spectrum become much more frag- 
mented. 



To get the two terminal conductance for a finite size 
diamond- Vicsek fractal, we attach the system between 
two semi-infinite one-dimensional single channel waveg- 
uides. The wave equation obeyed by the incident wave in 
these waveguides is discretized, and the leads are artifi- 
cially converted in to arrays of effective nodes character- 
ized by a constant on-site potential e; = 2 cos ka+2 cot ka 
as before, and a nearest neighbor overlap integral ti = 
1/sinfca. We then successively renormalize the finite 
network to reduce it into an effective two- vertex system, 
with renormalized effective on-site term equal to IA and 
with an effective hopping integral T between them. The 
transmission coefficient across the effective dimcr is given 



IV. CONCLUDING REMARKS 

In conclusion, we have examined the distribution of 
intensity of light (or any electromagnetic wave) on a Vic- 
sek geometry consisting of diamond shaped single chan- 
nel waveguides. The major result is that we have been 
able to identify a countable infinity of localized eigen- 
modes displaying a multitude of localization lengths. A 
prescription is given for an exact determination of the 
wave vectors corresponding to all such modes, a problem 
that is far from trivial in the case of a deterministically 
disordered system. The localized wave functions span the 
fractal space in clusters of increasing sizes, the size being 
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precisely controlled by the length scale at which the wave 
vector (k) is evaluated. The onset of localization can be 
exactly predicted from the stage of RSRG and can be de- 
layed (staggered) in space. The study provides a unique 
opportunity to experimentally observe the localization of 
light triggered by the lattice topology without bothering 
about the high dielectric constant materials. It may also 
be useful in developing novel photonic band gap struc- 
tures, where the wavelengths to be screened or allowed 
to go through can be controlled over arbitrarily small 



domains exploiting the fractal character of the network. 
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